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ABSTRACT 


A dynamic substructuring analysis, utilizing the component modes techni- 
que, of the 1/8 scale Space Shuttle Orbiter finite element model is presented. 
The analysis was accomplished in 3 phases, using NASTRAN RIGID FORMAT 3 (Level 
15.5.1), with appropriate Alters, on the IBM 360-370 (Model 165 ). The Orbiter 
was divided into 5 substructures, each of which was reduced to interface de- 
grees of freedom and generalized normal modes. The reduced substructures were 
then coupled in Phase 2 to yield the first 23 symmetric free-free orbiter modes. 
The eigenvectors in the original grid point degree of freedom lineup were then 
recovered in Phase 3. A comparison is then made with an analysis which was 
performed with the same model using the direct coordinate elimination approach 
under NASA contract NAS 1-10635-12 (Reference l). Eigenvalues were extracted 
using the inverse power method. 
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INTRODUCTION 


This portion of task NAS 1-10635-21 was undertaken to develop a modal 
synthesis approach to the sub structuring procedure for analyzing the 
elements of the NASTRAN finite element model previously generated for the 1/8- 
scale shuttle dynamic model. This model consists of an orbiter and two solid 
rocket boosters all attached to a central external tank. Photographs of the 
assembled model are shown in Figs. 1 and 2 (NASA Langley photos L73 6687 
and L73 6688). The NASTRAN (NASA Structural Analysis) finite element represen- 
tation of the orbiter model is described in Reference 1 . The NASTRAN finite 
element representation for the external tank and solid rocket boosters are des- 
cribed in References 2 and 3, respectively. A statistical description of 
these finite element models is shown on Table 1, which lists the number of grid 
points used, the number and types of members, and the degrees of freedom (DOF) 
remaining after reducing the number of independent coordinates. 

This reduction is accomplished by imposing single point constrants (SPC) or 
multiple point constrants (MPC), or by assuming certain coordinates have no 
forces applied to them. The latter approach is called Guyan, after its origi- 
nator (Reference k) . 

The overall analysis flow, in Fig. 3-1 in Volume II of Reference 1 , 
represents the originally proposed analysis for the combined total vehicle. 

The Orbiter was divided into five substructures: fuselage, cargo doors, fin, 

wing, and payload. The external tank was divided into two substructures: the 

LOX tank and the aft portion of the external tank (consisting of the inter- 
bank skirt, LH 2 tank, and aft tank skirt). The SRB originally was to be 
handled as a single unit (consisting of the forward skirt, propellant cylinder 
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and propellant, and the aft skirt), however, after computer storage problems 
were encountered, it was divided into two substructures as shown in Figs. l6 
and 17 of Reference 3. 


Referring to Fig. 3-2 in Volume II of Reference 1, observe that each of 


Llit; five Orbiter substi-uctui-es 


was analyzed to produce reduced mass and stiff- 


ness matrices for selected dynamic degrees of freedom (DOF's) and interface 
attachment points. Modes for these substructures were then obtained with the 
interfaces held. An exception is the fuselage, which was analyzed in a free- 
free condition. This approach aided in checking and understanding the behavior 
of the combined Orbiter vehicle. Next, the five substructure stiffness and mass 
matrices were merged to form the total Orbiter mass and stiffness matrices. 

These matrices were again reduced to yield final stiffness and mass matrices 
that were used in the modal analysis. This procedure of first merging the 
mass and stiffness matrices, then obtaining the eigenvalues, is called the 
direct method in this report. 


Several technical problems arose during the study which prevented the com- 
pletion of the proposed overall analysis, namely: 

• The Orbiter analysis was completed at the same time that initial 
test results were made available. A rather poor correlation was 
shown to exist for the Orbiter alone 

• The computer time required to analyze the hydroelastic model for 
the External Tank proved to be excessive 

• The computer time required to analyze the viscoelastic model for 
the Solid Rocket Booster as a single model was high. 
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Fig. 2 Assembled 1 /8-Scale Shuttle Model 
(Side View) 
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*In direct approach springs were included in coupling run. 































































































The two major problems encountered (lack of correlation of analysis and test 

data for the orbiter; excessive computer time requirements for cou pling the total 
vehicle) forced a decision to abandon the original overall analysis flow, con- 
sequently, baa|e. effort was redirected to rectifying the Orbiter analysis to 
obtain correlation with test results. The analytical and experimental investi- 
gation8 un ^|^i^te«cr I bed' in References 1 and 5. These resulted in re- 
vised orbiter f^aite ^lament representations which provided good agreement 


between a naly*s «* t*et. ^ response to the problem of excessive computer time 
a two-pronged study was undertaken under task NASl-10635-21 to find a means for 
improving the *^cien^y of the hydioelastic analysis and to develop procedures 
fbr using modal coupling for combining the NASTRAN substructure models. The 
latter effort isthe subject of this report# 

Much of th^teroi^lo®r describing the work done herein originates in the 
NASTRAN system, §p is picMhed in detail in Reference 6. 



ORBITER FINITE ELEMENT MODEL 


The Orb iter finite element model used in the analysis was the Model II ver- 
sion developed in NASA contract NAS1-10635-12 (Reference 1). The Orbit er was 
divided into five substructures (fuselage, wing, cargo doors, fin and payload). 
The Model J§ statistics on number of GRID points and types of finite elements 
are listed || Table 1. Also in Table 1 are the degree of freedom statistics 
for the modal synthesis and direct approaches. The NASTRAN Bulk Data for the 
various substructures are listed in Volume II. 




SUBSTRUCTURING PROCEDURE 

1116 **> the analysis of the Orbiter is 

known as the cotfegnt erodes or moded synthesis approach. The general theory 
is presented in Appendix A. The technique employed is essentially the same as 

w waxen is s imitar hut not identical to 

• ■' •;• "K; •'••;. 

Hurty s method. .• in . Ite?lireaoe •• B . The type of analysis chosen uses con** 

straint and noi*i^no4ep «w||tip»«3y, and the eigenvectors need not be normalized 
in any To a more reliable analysis, procedures 

to assess the validity of the steps in assembling the model were incorporated in 
the Direct Matrix Abstraction Procedure (DMAP) alter statements. These checks 
are inserted to insure that the oqnetraints applied by MPC’s and SPC s (multipedht 
and single point constraint#) d» not induce spurious loads or reactions into the 
structural ' aiod^^p9fti^ ; :^re to demonstrate that SPC’s do not' 

result in $1*8' o^i^ ; da#daw transformation matrices, such as g 

are checked to SflS&i^; th|pi is any deterioration in accuracy due to round-off 
or ill-conditioning . The reduced stiffhess matrix (after reduction) is checked 
for equilibrium. The reduced mass is converted to a rigid body mass (or weight) 
matrix so that it; eould be compared to the original matrix (MO matrix which is 
OUtpUt Generator)) before reduction. A more 

detailed description of the checks is presented in Appendix A. The NASTRAN steps 
for these pasoeeftii|^^fiE^'. 'itaom '■ ±n Appendix Bl. 

The theory was lncoiporated into NASTRAN Rigid format 3 via Alters. A de- 
tailed descriptiot'of the Alters can be found in Appendix Bl, while the actual 
IBM listing of Alters are in Volume II* Appendix B2. The analysis was performed 


in three phases, as shown in Fig. 3» for the schematic diagram of the analysis 
flow. The three phases are similar to those proposed by R. Guyan in Reference 9. 


A brief description of the three phases is as follows: 

• PHASE 1 - Component modes with interface fixed are calculated. The 

interface degrees of freedom are defined on SUPORT cards 
(r-Set). The interface supports can be determinate or in- 
determinate . Calculation of component generalized and 
reduced interface matrices (stiffness and mass) are per- 
formed and put on tape. Phase 1 is done for each Sub- 
structure. 

• PHASE 2 - In this phase all uncoupled interface points are defined 

on GRID cards. The same GRID cards from Hiase 1 can be 
degrees of freedom except at the interface 
are defined on SPC cards. All cooiponent inodes found in 

i' are defined by unique scalar point numbers* Higher 
^m : treqjuidcy inodes not considered essential can be put on SPC 
cards. The generalized and reduced interface matrices from 
Hiase 1 runs are then merged into an uncoupled pseudo- struc- 
ture -g lineup. The g-set consists of (6 x GRID POINTS + 

of freedom. The common Interface de- 
grees of freedom are coupled using MFC cards. The coupled 
structure can now proceed through the normal reduction pro- 
cess to yield system normal inodes. A tape is created for 
each substructure containing final eigenvectors in the sub- 
structure lineup which will be input to Hiase 3. 
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Fig. 3 Flow Diagram for NASTRAN 

Substructuring (Component Modes 
Method) to Obtain Orbiter Normal 
Modes 


















• PHASE 3 
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- Retrieval of final detailed substructure mode shape (eigen- 
vectors) in original substructure GRID POINT designation. 
Hiase 3 is done for each substructure. 
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RESULTS AND DISCUSSION 


Results of the analysis are presented and discussed in this section. Where 
possible, the results were compared to results obtained for the same model using 
the direct coordinate elimination approach, which was performed under RASA con- 
tract RAS 1-10635-12 (Reference I). 


PHASE I COMPONENT MODES RESULTS 

Although component modes were obtained in the direct method for checking 
purposes, the modes for this phase were either free-free, or restrained but 
without including the effects of the interface 'springs. Therefore, the only 
substructure that could be compared was the wing, and here the difference in 
frequencies was less than 1 %. 

Initially, in the modal approach, all but the massless degrees of freedom 
were retained to obtain component modes. The wing, which was the first sub- 
structure analyzed with this approach, yielded modes considered spurious. For 
example, Fig. k shows a fictitious mode caused by retaining dynamic coor- 
dinates at grid point directions connected by the minimal rods. The minimal 
rods were provided to prevent singularities in the idealization, since they 
separated shear panels not capable of resisting direct stress. These modes 
disappeared when the appropriate degrees of freedom were omitted by GUYAN re- 
duction. 

Omitting only the massless and fictitious degrees of freedom worked well 
with all substructures, except the fuselage. Here, additional coordinates had 
to be omitted, since the number retained in the direct method was an upper 
limit, if excessive computer time was to be avoided. 
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MODAL OtTOR. SUtCASC 4 MOOC 4 FRCQ. 2CO.3S30 

- ' ’ • ' 1 ! 


Fig. 4 Fictitious Wing Mode Caused by Not 
Omitting Degrees of Freedom in 
Direction of Minimal Rod Line 
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Component mode plots are presented in Volume II. They contain 57 fuselage 
modes, 20 wing modes, 35 cargo door modes, 7 fin modes and 12 payload modes. 

A closer examination of these plots uncovered some deficiencies in the Model II 
idealization of the wing. 

The 6th wing component mode (404.5 Hz) in Appendix B6 demonstrated that 
this idealization had practically no lateral resistance at the interstage 
station. These flaws are shown in Fig. 5, which also indicates the fix-up 
to be taken. This error should not affect the total Orbiter system modes, but 
it would certainly affect a total Shuttle analysis, where the inclined inter- 
stage link would produce force componenets in the lateral direction. 

The 10th and 11th wing component modes (599. 4 and 613.6 Hz) in Appendix 
B6 revealed the other flaw indicated in Fig. 5* The above modes disappeared 
when the wing was rerun through Riase 1 with the indicated modifications. Table 

2 contains comparison of frequencies before and after the modification. Figures 

6 through 15 show plots of the modes of the revised wing. Comparison of the 
modal plots before and after modification show that the "kinks" have disappeared. 
The revised wing was not used in the Orbiter analysis, since comparison of re- 
sults with the same model that was used in direct method analysis was the objec- 
tive. The final orbiter results (first 23 modes) indicated that only the first 

3 wing component modes played a significant part for most orbiter modes (Refer to 

Table 5). The small difference (3$) in frequencies for the first 3 modes 

(Table 2) would not have influenced the Orbiter results appreciably. 


SOLID LINES AEE GRID LINES WITH STRUCTURAL RODS 

DOTTED LINES ARE GRID LINES WITH MINIMAL RODS (A = .001 IN - ) 
SHEAR ONLY PANELS 


FIXUP (BOTTOM COVER ONLY) 

(?) SUB-DIVIDE PANEL & INSERT 
STRUCTURAL ROD BETWEEN 


NO FINITE ELEMENTS 
IN TRIANGULAR 

REGIONS (TOP & BOTTOM COVERS). 

ANY SHEAR IN PANEL NEXT 
TO MINIMAL ROD AND TRIANGULAR 
HOLE, WILL DUMP LOAD INTO 
MINIMAL ROD WHERE THERE IS 
PRACTICALLY NO RESISTANCE 

FIXUP (TOP & BOTTOM COVERS) 
(3) INSERT STRUCTURAL RODS 


PRACTICALLY NO 
RESISTANCE TO 
FORCE IN DIRECTIO 
SHOWN 



Fig. 5 Flaws in Model II Idealization (Bottom Wing Cover Shown) 
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Fig. 6 Revised Wing (Mode 1 
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PHAMC 1 (OftfttTCft W1N®> 

via n« <cavc** to pcrccnt crr.> 
mcc mooc« fixed at interface 
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Fig. T Revised Wing (Mode 2) 
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Fig. 8 Revised Wing (Mode 3) 


PHABC I CORBITT ft WIM9> 
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Fig. 9 Revised Wing (Mode k) 
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Fig. 10 Revised Wing (Mode 5) 
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Fig. 11 Revised Wing (Mode 6) 
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Fig. 12 Revised Wing (Mode 7 ) 
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Fig. 13 Revised Wing (Mode 8) 
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Fig. Ik Revised Wing (Mode 9) 
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Fig. 15 Revised Wing (Mode 10 ) 


Table 2 Wing Substructure Component Modes Comparison 
of Model II (Before and After Fix-up) 




FREQ* 

AFTER 

FIX-UP 

HZ 

COMPONENT MODE DESCRIPTION 

1 

74. 5 

74.8 

1 S ^ Vert. Bend. 

2 

148.5 

153.3 

_ st _ 

1 Torsion 

3 

254.9 

259.7 

_ st 

1 Lateral Bend. 

u 

298.4 

304.0 

2 nd Vert. Bend. 

5 

330.4 

340.8 

0 nd m 

2 Torsion 

6 

4o4.5 

- 

local Lateral at Interstage due to Idealization Error 

7 

517.0 

526.9 

Combined 2 n ^ Vert. Bend & Torsion 

8 

542.1 

564.1 

2 nc * Lateral Bend. 

9 


584.4 

Combined 2 nc * Lateral Bend & Vert. Bend 

10 

599-4 

- 

Caused by Idealization Error at Triangular Cutouts 

11 

613.6 

- 

Caused by Idealization Error at Triangular Cutouts 

12 

648.2 

652.9 

rd 

3 Vert. Bend. 

13 


679.8 

rd 

3 Lateral Bend. 


*Fix-up version not used in overall analysis 
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FINAL SYSTEM ORBITER RESULTS (SYMMETRIC MODES) 


Table 3 shows good correlation between the Modal Synthesis and direct 
method, and verifies the analytical eigenvalues for the 1/8- scale Model II 
Orbiter. The door modes which are apparent in the Modal approach were completely 
overlooked in the direct method because no dynamic degrees of freedom were in- 


uxuucu vxxi oxic: uvux i xixoxxuugxj the dOOx xuaSS 15 small, the duux' LOxigex'vil, which 
is the only door structure working in the symmetric case was flexible enough in 
bending to cause lower frequency modes. Table 4 shows the contribution to the 
Orbiter System generalized stiffness and mass of each substructure. As can be 
seen, the door contribution is small for the overall Orbiter modes. Table 5 
shows the contribution factors (generalized modal coordinate values) of the 
substructure component modes to the Orbiter Sym. free-free modes. Table 6 
gives the major contributing component modes and an attempt to classify them. 
Plots of the final Orbiter modes were obtained from a Phase 3 run for each sub- 
structure. Therefore, plots of the total Orbiter on one sheet could not be 
obtained. Instead, individual substructure plots were obtained which were of 


different scale. The plots are shown in Appendixes Bl4 through Bl8. As a 
final check. Table 7 shows the summation of relative momentum forces of each 
substructure for each mode and should demonstrate momentum balance (Sum = 0 
for the Orbiter). 


Because the door motion was prominent in the combined Orbiter modes, it 
was decided to rectify the discontinuity between the door longeron and shell 
at the interior strap locations. This can also be seen in the Model II door 
component mode plots (Appendix B8). This was done by reverting back to a 
Model I door shell where a shell grid line was incorporated at the interior 
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strap stations. The subdivided panels were made of CQUAD2 elements. This re- 
vised door was run through Phase I to obtain component modes. The results of 
this mm are presented in Table 8 and the modes plotted in Figs. 1 6 through 
27- From Table 8 the sensitivity of the door can be seen, when some of the 
skin is made effective to resist the door longeron bending. The frequencies on 
the average were higher and the sequence of some of the modes was changed. This 
revised door was not incorporated into the overall analysis since, as previously 
stated, the purpose was to compare two methods for the same model. 

Actually, the door should be idealized into a finer grid (station-wise) to 
correctly represent the skin contribution to the bending stiffness of the door 

longeron, which will materially affect the mode frequency. This would also help 
if an Orbiter anti - symmetric analysis were to be performed, where the shear in 
the door panels is transferred through the door longeron and straps to the 
fuselage. The anti-symmetric torsion test case (Reference l) on the l/8-scale 
Orbiter proved that the Model II door longeron was too flexible since the analy- 
sis gave twice the torsional influence coefficient. 
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Table 3 Comparison of Analytical Results Between Substructuring Methods 
for Symmetrical Free-Free Normal Modes (l/8-Scale Model II ) 
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Table U Substructure Contribution to Generalized Stiffness and Mass 
of Orbiter for Symmetric Free-Free Modes (l/8-Scale Model 
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substructure generalized stiffness matrix = 0 t k0 
substructure generalized mass matrix = 0 T M0 




















Table 5 Contribution Factors (Generalized Modal Coordinate Values) of Substructure 
Component Modes to Orbiter Symmetrical Free-Free Modes (l/8-Scale Model II ^ 
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Orbiter Modes 1 to 3 are rigid body modes. 




























































































Table 6 Substructure Component Modes (Symmetrical 
Case) l/8-Scale Model II 


SUBSTRUCT. 

MODE 

NO 

FREQ. 

HZ 

MODE 

DESCRIPTION (FIXED AT INTERFACE) 

Fi 

iselage 

1 

86.9 

Nose Pitching Mode 






2 

162.5 

Nose Fore -Aft Translation Mode 




3 

245.6 

Mid Section 1st Bending Component 

Mode 



4 

270.6 

Aft End Vert. Bend. 






5 

280.8 

Mid Section 2nd Bending Component 

Mode 



6 

333.8 

Nose Vert. Bend. 






7 

339-9 

Aft Frame Mode 







8 

378.9 

Mid Section 3rd Bending Component 

Mode 



9 

391.6 

Mid Section 4th Bending Component 

Mode 



10 

439.8 

Aft End Vertical Translation Mode j 



11 

448.2 

Mid Section 5th Bending Component 

Mode 



12 

498.0 

Local Thrust Bar Translation Mode 


f 

17 

633-8 

Nose & Mid Section Bending Mode 




24 

839.0 

Local Thrust Bar Axial Mode 


Wing 

1 

74.5 

1st Vert. Bending 






2 

148.5 

1st Torsion 







3 

4 

254.9 

296.4 

1st Lateral Bend. 
2nd Vert. Bending 


Cantilevered Modes 



| 5 

330.4 

2nd Torsion 







8 

542.1 

' 2nd Lateral Bend. 




Cargo Drs. 


47.5 

1st Vert. Bend. 








53.5 

2nd Vert. Bend. 



Continuous Door Longeron Beam 




79.7 

3rd Vert. Bend. 



Modes on 

6 Supports 




89.9 

4th Vert . Bend 







5 

130.5 

Local Aft Frame Mode 





6 

131.0 

Local Fwd Frame Mode 





7 

163.6 

1st Lateral Bend. 






8 

174.1 

2nd Lateral Bend. 


Continuous Door Longeron Beam 



9 

190.8 

3rd Lateral Bend. 


Modes on 

Flexible Frame Supports 



10 

226.0 

4th Lateral Bend. 




\ 


11 

314.9 

Fore -Aft Translation 

Mode (Supported at strap interface) 



12 

477.7 

Fore -Aft Axial 

Mode 




Fin 

1 

77.6 

Pitching Mode 

NOTE: Fin on vert, interface springs 

1 

2 

420.9 

Vertical Mode 





Payload 

1 

647 ~ 

1st Vert. Bending 

Vert. Bend. ) 

Simple supported modes 



2 

131.2 

Pitching (some 

1 st 

on vertical interface 

■ 


3 

163.3 

Pitching & 2nd 

Vert. 

Bend. 

1 springs 

i " 

4 

373.0 

2nd Vert . Bending 
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Table 7 Summation of Substructure Momentum Forces About Basic Origin 
For Orbiter Symmetric Free-Free Modes (l/8-Scale Model II) 




















Table 8 Cargo Door Substructure Component Modes (Symmetrical 

Case) Comparison of Model II (Before and After Fix-up) 


DTP TTirvnTT* 
XJXJ J.' UillJ 

T7>TV TTPi 

r 1A- U JT 

AFj.uR 

FIX- UP 

DESCRIPTION 

MODE 

FREQ. 

HZ 

MODE 

FREQ. 

HZ 

1 

47.5 

i 

50.2 

st 

1 Vert. Bend Component Mode 

2 

53.5 

2 

58.5 

2 n< ^ Vert. Bend Component Mode 

3 

79.7 

4 

124.1 

rd 

3 Vert. Bend Component Mode 

4 

89-9 

3 

116.1 

t h 

4 Vert. Bend Component Mode 

5 

130.5 

5 

131.0 

Local Aft Frame Mode 

6 

131.0 

6 

131.2 

Local Fwd Frame Mode 

7 

163.6 

10 

388.9 

s t 

1 Lat. Door Long. Bend. 

8 

174.1 

n 

450.2 

2 nC ^ Lat. Door Long. Bend. 

9 

190.8 

7 

200.4 

rd 

3 Lat. Door Long. Bend. 

10 

226.0 

8 

228.0 

4 th Lat. Door Long. Bend. 

11 

314.9 

9 

316.8 

Fore- Aft Tran s la t ion 

12 

477.7 

12 

487.6 

Fore-Aft Axial j 


*.Jote: Fix-up version not used in overall analysis. 
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Fig. 1.6 Revised Cargo Door (Mode l) 
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Fig. 17 Revised Cargo Door (Mode 2) 
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Fig. 18 Revised Cargo Door (Mode 3) 
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Fig. 19 Revised Cargo Door (Mode 4 ) 
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Fig. 21 Revised Cargo Door (Mode 6) 



Ui 


Fig. 22 Revised Cargo Door (Mode 7) 
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. Fig. 23 Revised Cargo Door (Mode 8) 



>4 3 
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Fig. 25 Revised Cargo Door (Mode 10 ) 



^5 


i 

i 


Fig. 26 Revised Cargo Door (Mode 11 ) 



. 27 Revised Cargo Door (Mode 12) 


COMPUTING TIME 


It seems logical in the Modal Synthesis method to keep as many dynamic de- 
grees of freedom as possible in obtaining component modes. Thus approximations 
brought about by GUY AN reduction are eliminated. This approach worked rather 
well with the smaller substructures, where only massless degrees of freedom and 
those along fictitious nodal lines were eliminated by GUYAN reduction. The 
large fuselage substructure could not be treated in the same way. It took 24 
Central Processing Unit (CPU) minutes or 284 system minutes to obtain 1 mode with 
6l4 degrees of freedom. Nine (9) (CPU) minutes or 115 system minutes were spent 
in the Real Eigenvalue Analysis (READ) module alone. It was therefore decided 
to keep approximately the same number of dynamic degrees of freedom that was used 
in the direct approach to compute component modes. Although there was no choice, 
this was compounding the lack of accuracy since in modal synthesis accuracy is 
lost by carrying only a reduced number of component modes into the coupling 
run. Table 9 shows a comparison of computing time to obtain Orbiter Symmetric 
modes for the two methods. The time spent in the READ module in system minutes 
per mode extraction as a function of the dynamic degrees of freedom is plotted 
in Fig. 28. As can be seen in the figure, there is a great need to incorporate 
into NASTRAN a more efficient eigenvalue extraction program, especially if one 
has to calculate higher modes for a large problem. For higher modes there should 
be more dynamic degrees of freedom (less Guyan reduction) in calculating component 
modes, and more component modes must be extracted in Phase 1. 

From Table 9, the direct method is more economical, if one is to solve for 
only the lower modes of the Orbiter. The advantage of the component modes ap- 
proach lies in the Phase 2 or coupling run, if more substructures were to be 
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coupled to the Orbiter (total Shuttle), For example, the Phase 2 results (Refer 


to Table 5) shows that only about ^0 component modes were important in computing 
the first 23 system modes. The unimportant component modes in Phase 2 could be 
eliminated, thus reducing the dynamic degrees of freedom from 220 to 121. From 
Fig. 28, it is evident that eigenvalue solutions in relatively short time can be 
obtained up o 0 about tbO degrees of freedom. This leaves approximately 250 de- 
grees of freedom for the added reduced substructures (external tanks and SRB) to 
be solved within a reasonable time for the real lower modes. 



Table 9 Computing Time to Obtain Orbiter Symmetric Modes 






OBSERVATIONS AND RECOMMENDATIONS 


The three phase modal coupling procedures adapted to NASTRAN may be sum- 
marized as follows: 


• The finite element model is divided into convenient substructures. 


ATT 4 -^ 4 - « J — — _ -T» 3 1 l — ~ T — 1 1 — -? n • — 

±uociiav>c uc^icco iiccuum ucuwcca cj-Giueu oc> clj. c it: uctJLiitfU. J_ll 


the analysis 


• Phase 1 consists of calculating modes for each substructure re- 
strained at the interface points. This phase also includes a series 
of check calculations to guard against spurious forces, loss of mass, 
or ill- conditioning during the matrix reduction process 


• Phase 2 consists of assembling the substructure models together. In 
this phase the interface points are merged and the modes considered 
unnecessary to represent the various substructures are eliminated 

Check calculations are made to test the merged interface stiff- 
ness and mass matrices. The eigenvectors and eigenvalues for 
the combined system are calculated 

• Phase 3 consists of retrieval of the final detailed mode shape 

for each substructure from individual tapes prepared during Phase 2. 

This method was applied successfully to the l/8-scale shuttle model orbiter. 
Comparison between results from modal coupling and those from the direct sub- 
structure merging method previously used indicated good agreement. Spurious 
modes of the cargo bay doors found here were probably also present in the 

other method but were not revealed since no dynamic degrees of freedom were 

assigned to them. 
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The method can be used to reduce the dynamic degrees of freedom for the or- 
biter, and to add the modes of the external tank and solid rocket booster while 
limiting the problem size to about 350 degrees of freedom, which should permit 
a real eigenvalue analysis of the combined shuttle. 

Using Modal Synthesis and a harmonic reduction technique developed by Robert 
Coppolino (Reference 2) for the hydro-elastic tanks, it is possible to reduce 
the final dynamic degrees of freedom for that substructure down to approximately 
350. Lower real eigenvalues could then be extracted within a reasonable time. 

Regrettably, the total l/8- scale model of the shuttle could not be analyzed 
because of limited time available, and only the Orbiter was analyzed to test 
the Modal Synthesis procedure developed in this report. 

Damping was not included in the modal synthesis procedure, since it is 
only considered significantly large in the SRB substructure, which was not in- 
cluded in this analysis. This would necessitate incorporating modal synthesis 
procedures into Rigid Format 7, which uses the complex eigenvalue module CEAD. 

It is doubtful if the total Shuttle model could be reduced to only 150 meaningful 
dynamic degrees of freedom, as required, to have the complex eigenvalues ex- 
tracted within a reasonable time. The 150 figure is based upon experience in 
analyzing the Solid Rocket Booster ( SRB ) , Reference 3. Twelve (12) complex eigen- 
values for 116 dynamic degrees of freedom were obtained, which took about 6.3 
system minutes per mode in the CEAD (Complex Eigenvalue Analysis) module . This 
is about six times as long as a comparable problem required in the READ (Real 
Eigenvalue Analysis) module. Referring to Fig. 28, we can see the system 
minutes per mode for the CEAD module will rise rather sharply compared to the 
READ module as the number of dynamic degrees of freedom increases. 
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PER MORE 


It is therefore strongly recommended that before a modal synthesis (sub- 
structuring) procedure is adopted to yield complex eigenvalues, a more efficient 
complex eigenvalue extraction program be developed. 



Average -Lue f'pent in READ Module Extracting 1 Mode 
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Appendix A 

NASTRAN COMPONENT MODES ANALYSIS GENERAL THEORY 


APPENDIX A 


r 


APPENDIX A NASTRAN COMPONENT MODES ANALYSIS - GENERAL THEORY 

Phase 1 - Representing Part of a Substructure by Normal Vibration Mo des 

The equations of motion for a substructure (after GUYAN reduction, if any) 

are 

M ■ |V] N + [ K aa] {“a} l) 



Or, letting jp^} represent the vector of Inertia forces, then 


[ K aa] j Ua | ' j P aj 2) 

where 

H " H |“a| 3) 

Bartition eq. 2 into interior and interface degrees of freedom ( land rsets) 



K 


1 - 


£r 


K J 
rr 



4 ) 


The substructure displacements may be represented as the superposition of 
displacements relative to the interface and those due to interface motion, as 
follows : 


A-l 


5 ) 


i u 4 = W> + 

o s, > 


where ^u^ is the vector of displacements relative to u r (i.e., with = ^), 

and ^u^ is the vector of displacements due to ^ u r^ . The ^_ u j^ displacements are 
due to £pj?\ with ^ u rl, = l_ 0 } , while the ^u^ displacements are due to ^ u r "^, with 


$pJl = $ o\ . 

c ^ ^ 


The relationship between £u £ ^ and ^u r ^ may be determined from the upper parti- 
tion of Eq. 1+ as 


& - M K'5 


6) 


where 


l G Jt>'L K u} X W] 


7) 


combining eq. 5 and 6 gives 


_ u £ 

u„ 


I _^£ 

! Ir 


U£ 


U T 


8) 


where I and I are unit matrices. 

£ r 

Partitioning eq 1 into interior and interface degree of freedom, gives 


ffj 'v vv ; 

I \= CVT * v+ 


'• K ! K i U 

i ££ Jlr £ 


9) 


... V l t 
1. . I v K 


(r rl ; v i VI C r ) 


U 


rr 


where 


'h > 

0 


A-2 


holding the interface fixed and writing the upper partition of eq. 


relative 


displacements 




, gives 


9 for the 


+ [k£jQ 10 ) 

The corresponding real eigenvalue problem is 



letting 


to -wto 

where 

= [ i 2 — [ 0 ^ b~\ 


and 


&v 


vector of modal displacements, Eq. 8 may be rewritten as 


12 ) 



13) 


It should be noted that^0jQ contains a reduced number of modes, i.e., the 
number of columns of to is less than £ , which is the number of degrees of 
freedom in eq. 10. 

The generalized modal forces, as shown on page 14.1-3 of Reference 10, can 
be expressed as: 



A-3 



Substitution of first Eq. ik then Eq. 13 into Eq. 9, and using Eq. 7, yields 
the following reduced matrix equation of motion in terms of the generalized modal 
and interface coordinates. 




15) 


where 



- W'WW 

16 ) 

l K rr") 

= [%rlV<f] + t rr l 

17) 

[mu] 

-M [M U ][^ 

18) 

[%rl 


19) 


+ £ M i.r]jj + [%r T ][pJ2.1 + [_Mrr~] 20) 


The resulting matrices from Eq. 16 to 20 will be input to Phase 2, to be 
coupled to other substructures. 

Kn and are the generalized modal stiffness and mass matrices and 
should be diagonal matrices. In the actual NASTRAN computations, small off- 
diagonal non- zero terms will occur. The following relationship should exist 
between the corresponding diagonal elements 


k ii ~ ^ i m ii 


21 ) 
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K and M are the static reduced interface stiffness and mass matrices, 
rr rr _ 


when the interior degrees 


of freedom are released. would 


repre- 


sent the generalized modal forces, due to interface accelerations. 

PHASE 1 - INCORPORATED CHECKS 

The checks incorporated are all based on a matrix EQg, which can be ex- 
tracted from the GIWG module. The NAS TRAN module GFWG has been modified (Appen- 
dix B2) to output this matrix, which expresses the static load summations for 
each unit g-set load. This procedure is forced when the parameter WTMASS, in 
the general calling sequence of GIWG (NASTRAN Programmer's Manual (NPM) 4.29), 
is set to 0.0. For all other values of WTMASS, module GIWG performs as outlined 
in the NH4. EQg is a 6 x g matrix, where g is equal to 6 times the number of 
grid points in the problem. It should be noted, that only grid points should be 
used in the problem when extracting this matrix, since scalar points have no 
geometry. Therefore, in Blase 1, only grid points are used. The 6 rows of EQg 

correspond to the 2F , 2F , 2F , 2M , 2M , and EM load su mm ations respectively 
* x’ y’ z’ x’ y z 

about a reference point specified by the parameter GRDPNT. An example of ex- 
tracting EQg from GIWG is as follows: 


GPWG BGPDT , CSTM , EQEXIN , /EQg/ V , Y , GRDPNT— l/C , N , 0 . 0 


where 


i F REF^ ■ EQ g i 

__6xg 


STATIC EQUILIBRIUM 


As indicated for equations 13 and l4, the following transformation holds 


KINEMATIC CONTINUITY 



where 


M ■ UJ 


24 ) 


and 


VUrEf"^ is the vector of 6 rigid body motions of the reference point and 
^U^contains all g-set displacements. 


Matrix D 
Manual , 


is equal to the D matrix discussed in the NASTRAN Programmers 
g 

Section 4 . 29 . 


Matrix D can he partitioned into the various NASTRAN subsets by using 

g 

column partitioning vectors generated by the VEC instruction. The subsets of 

D are as follows 

g 



MULTIPOINT CONSTRAINT CHECK 

The NASTRAN program forms the matrix Gm from the MPC bulk input. 


^ u m Gm \ ^ u n^ 


29 ) 


A-6 


can be related to rigid body motion at the reference 


The displacements 

point by D in Eq. 25. 
n 




“ D n~] £ u REf£ 


30) 


Equation 30 should be equal to 


£ u m"^ "[ D m 1 ^ u REf} 


or 

(toco - cbj)[u r 4 = H 


31) 


or 


[MPCCK] =[G m l[ D J = [o] 

mx6 


32) 


When performing symmetrical or anti -symmetrical analyses, MPCCK may contain 
non-zero terms if the dependent degree of freedom is on the plane of sy mm etry 
and the independent degree of freedom is off the plane of symmetry. In this case, 
the non- zero term will be a difference in coordinates between the 2 points. 

For example: 

PLANE OF SYMMETRY' 


u_ ( INDEPENDENT ) 
2z 

t 

u— T — 

z 


. u (DEPENDENT) 

f 



REF. 


PT. 


I 
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9 


u iz 


u 

zz 


u. 

lz 


= t G m *1 Ufcz 

lxl 



[mPCCk! =lG m l[D n l 


CO - Cil 

Hr) n l = U 0 I o' l| Y2 i o! o'] 
C = L° 1 ° \ 1 1 ° i o l o*] 
- LdI =[oi o; 0! Y2| oj o'] 

*“ m J * 

^ Non-zero 


The six columns pertain to rigid "body motion at the reference point 

( u ox , "oy > > °ox > ©oy > ©oz) • The non-zero term is caused by an anti -symmetric 
motion ©ox, which doesn't apply to the illustrative symmetrical case. The terms 
under the symmetrical motions Uo X , u oz , & 0 O y in this case should always be 
zero . 

The succeeding checks to be developed on the following pages will also follow 
the same rule. That is, when performing symmetric or anti- symmetrical analyses 
only the three related columns are appropriately looked at in the check matrix. 

SINGLE-POINT CONSTRAINT (SPC) CHECK 

An SPC check is developed which is based upon the following assumption. The 
only degrees of freedom to be included in this set will be those that have no 
stiffness and those that are symmetrical or anti -symmetrical boundary constraints 
at the plane of symmetry. Any other supports that a structure might have are 
included in the r-set (SUFORT card). Appropriate ALTERS, to change the condition 
that the r-set by statically determinate, have been made to prevent a FATAL ERROR. 
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The following matrix is formed in NAS TRAN 


w =ro T h 


33 ) 


The displacements ^_ u f"^ can be related to rigid body motion at the 
point by[D f "]in eq. 26. 


reierenee 


to - L spcck 3 Kef^ 


34 ) 


where 



35) 


QSFCCK should be null. For symmetrical or anti symmetrical analyses only 

the appropriate three columns will be zero. 


SINGLE- POINT CONSTRAINT MASS CHECK 


When mass is generated from member densities, mass may inadvertantly be 
assigned to SPC degrees of freedom. This mass will be lost in calculating vibra- 
tion modes, unless they happen to be at sym. or anti boundary constraints. 

SPC inertia forces can be written as: 


to - Kslkl 36) 

where £ M SS ~3 is a symmetrical partition of £ M nn l * 

The accelerations can be related to rigid body accelerations at the 

reference point by [_D S ~] in eq. 26. 
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{f s \ = QmspcI ^ref} 


37) 


where 

EmSFc! = H“ssl [ D s] 38) 

5x6 

U MSPC should be null. For symmetric or anti -symmetric analyses only the 
appropriate 3 columns will be zero. If they are not, the degree of freedom in 
question should be MPC'ed, to prevent loss of mass. 

OTHER TRANSFORMATION CHECKS 

Checks similar to the MPC check (eq. 32) are performed for the NASTRAN 
generated transformation matrices C G 0 ~] and£(j , where [ G^= . 

This was done mainly to determine how far equilibrium has deteriorated due to 
ill-conditioning or round-off. The checks are: 


QO 1 = [pOCHK~] = [gJd^ - [V] 

39) 

[ 0~] =[GLCHK"\ = [G ? ][D r l - [d £ 1 

40) 


For symmetric or anti -symmetric analyses only the appropriate 3 columns 
should be zero. 

REDUCED INTERFACE STIFFNESS CHECK 

The static interface stiffness from eq. 17 states 




4l) 
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Relating £u r l- to rigid body motion by ^_ D r3 (eq. 28). 

t°l - &rrck] I* REr'i =[K rr lLD r l **2) 

^ Null 

For symmetric or anti- symmetric analyses only the appropriate 3 columns should 
be zero. 

RIGID BODY MASS MATRIX CHECK 

The reduced interface mass can be converted to a rigid body mass matrix. 
This can be compared with the £mcQ matrix, which is printed output from the 
GFWG module. For symmetric or anti- symmetric analyses, only the symmetric or 
anti-symmetric terms should be compared. This check ensures that no mass has 
been lost in the reduction process, is converted to a rigid body matrix 

as follows: 


QmokrI = 01 T , l M 1 » 01 43, 

6x6 ix 

MATRICES GENERATED IN PHASE 1 NECESSARY FOR PHASE 2 CHECKS 

In Phase 2, the basic matrix EQg (eq. 22) cannot be extracted from the 
GfWG module, because the Riase 1 component modes (or generalized coordinates) 
will he defined in Riase 2 as scalar points. Therefore, it is necessary to 
generate matrices in Riase 1 which can be used for Riase 2 checks. 

We already have a matrix [iT](eq. 28) to define the interface motion due 
to rigid body motion at the reference point. This matrix will be input to 
Riase 2. 
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kb) 


i u r*^ ~ C Dr "1 £ u REF^ 


We must now find a similar matrix for the generalized modal coordinates, 
which will be written as 


= [DilfupiiJi 

' ✓ V -/ 

Inverting \Kii^ in eq. 15 yields 

N - { p t ] 

The generalized forces, ^Fi"^ , is defined in eq. lil- 


ts 


k6) 


as 


t^Vixy 1 w 

The generated NASTRAN matrix [k^"] defines 

[ F ^=[sulK\ m 

f • _ 

Converting to rigid body motion by^"} in eq. 28 gives 

= 1 U REF^ ^9) 

Sincejjfjj^ is large ^ F^ can be defined another way by using the upper 
partition of the stiffness matrix in eq. 9) 

i F ^ =l K ul + l K lr~} 50) 
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51 ) 


Converting the displacements to rigid body displacements will set ^F^ = ^0^ 


or 



LVlOll - - [%rH>r] 

Therefore, substituting 52 into 49 yields 


" t K irlt D r1 ^_ U REF^ 

Combining eq. 46, 47 and 53 gives 






Equating 54 to 45 yields 



52) 


53) 


54) 


M = - iHif [ K 1 l *irl£>r 1 55 ) 

This matrix will be input to Phase 2. 

The column partition vectors used for merging substructures in Haase 2 now 
seems to be the only unchecked hand data. These vectors are somewhat incon- 
venient to prepare and are subject to human errors. Certain matrices will now 

be generated in Phase 1 so that they can be compared with the merged matrices 
in Phase 2 . They are : 


[ F REF^ = [ D i3 T c K iil = £ SDMKiil ^ 

6xi 


56) 
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57 ) 


£ F REF^ = t D il T [ M ii] (? = C SUMMiil ^ 

M = ^ D rT [ M irf ^ = C SUMMri] 58) 


[sUMKii“] gives the summation of interior elastic forces about a reference 
point due to unit generalized modal displacements. 

^SUMMii"] gives the summation of negative interior inertia forces about a 
reference point due to unit generalized modal accelerations. 

|T SUMMri^ gives the summation of negative interface inertia forces about a 

reference point due to unit generalized modal accelerations. 


Phase 2 Coupli ng Substructures' Reduced Dynamic Equations and Solving for 
Free-Free Modes 

The equations of motion of the combined uncoupled substructures can now be 
written in the following form: 

C M ggl [ u gl + C K ggl f u g} = 59 ) 

where £ F g^ = 


or 


j-MGG 

r 

JGG^ 


MSGpi 

MGGi 



where 



0 

0 


1 



6o) 
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*i u \ represents all of the uncoupled interface or junction point degrees 
of freedom. The interface or junction points are defined by GRID cards, thereby 
creating 6 degrees of freedom at each junction point. The unwanted d.o.f. (those 
that have been eliminated in Phase l) will be SPC'ed out subsequently in the re- 
duct ion process. The lineup of w is 


/ i \ 

u- 


H - « 


where N = no. of substructures 


u 


,N 


would contains 6 x (number of substructure i junction pts). 

represents all of the uncoupled generalized modal coordinates ob- 
tained in Hiase 1. These coordinates will be defined by scalar points. Any 
unwanted generalized coordinate (those representing higher modes) can be SPC'ed 
out subsequently in the reduction process. The lineup of £^\is 


S 1 


isV 


N’ 


The column partition vectors used to merge the substructures can be thought 
of as transformation matrices for the sake of presentation. For example: 

W3- W3W1W 

*fch 

will merge the j substructure interface stiffness from Haase 1 into the Hiase 2 
g-lineup. Therefore, the mass and stiffness matrices in eq. 59 are generated by 

][T 6 n T + Ogi ] C M iiI T gi] T + [ fgllMir 

- 
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C K gJ' (L TgJ r][^][^f 

where N = number of substructures. 

Similarly the kinematic matrix similar to eq. 23) can be generated 

l D g] = ^1 (C T gr][4]*[Tgil[ 4]') 61) 

In order to partition the merged matrices of eq 59) into that of eq 6o), we 
need a column partition vector defining the generalized coordinates in terms of 
the g-set lineup. This is obtained by adding up the substructures’ partition 
vectors which merged the substructure generalized coordinates into the Phase 2 
lineup. 



We can now perform some checks on the merged matrices in eq. 60) a fter 
first partitioning 



MERGED INTERFACE STIFFNESS CHECK 

Co} = \_KGGRCk 1 = L KGG rl[ DG rl 64) 

For symmetric or anti- symmetric analyses only the 3 appropriate columns should 
be zero. 
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MERGED RIGID BODY MASS MATRIX CHECK 


T 

j^MOGG r ~] = [DG r "J £ M3G r 1C D G r l 65) 

this matrix should be equal to 


^mogg/] 



Q MORR^ "] 


66) 


or equal to the sum of the substructures' rigid body mass matrices given in 
eq. 43. 

FORCE SUMMATION CHECKS OH MERGED MATRICES 


[sUMKGGil = 
6x 5 

CsUMMJGi'] = 
6x § 

HsUMM3 ri "3 = 

§ 


Oil C. H»i 1 

L“ifC 

t re rf [. 



Merged Matrices 


r SUMKj_i j 

— I SUMkJJ *1 

fsUMMii 1 

1 N -1 

— \ SUMMii ^ 

L sDiw ri ; — ■ 

N 

— » SUMMpi J 


Phase 1 Matrices 
eq 56 — » 58 


PHASE 2 CONTINUATI ON 


After the merging checks are performed, the dynamic problem stated in eq. 

59) will be reduced in the normal RIGID FORMAT 3 fashion to obtain a real eigen- 
value solution. Continuity at the interface between structures are described by 
MPC’s and the rigid body supports described by a SUFORT card. The checks incor- 
porated in Phase 1 are incorporated in Phase 2 (Equations 25 thru 43 are still valid 
in Phase 2). 
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The system eigenvalues and eigenvectors are recovered in the original sub. 
structure lineups and put on individual substructure tapes so that grid point 
displacement can be obtained and plotted for the system modes in H,ase 3. 

The eigenvectors for a typical reduced substructure would be 



67 ) 


Each substructures' system modal stiffness and mass is also calculated and 

printed out in this phase. This gives us the contribution of each substructure 
to the total modal stiffhess and mass. 

For more detailed description of all operations performed in all three 
phases, see Appendix Bl. 
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Appendix B1 

NASTRAN COMPONENT MODES ANALYSIS ALTERS TO RIGID 
FORMAT 3- PHASES 1, 2, & 3 


APPENDIX B1 



APPENDIX B1 NASTRAN COMPONENT MODES ANALYSIS - ALTERS TO RIGID 
FORMAT 3, PHASES 1, 2, AND 3 

REGULAR BULK PARAMETER USED - PHASE 1 

GRDFNT This parameter should always be used. It causes the rigid body 

mass matrix MO to be printed out, which can be compared with 
the matrix MORR discussed in ALTER 75,84. 

WlMASS - - Converts generated weight to mass. In the 1/8 scale model, the 
weight was in lbs., therefore WTMASS = .002588. The MO matrix 

was thus a rigid body weight matrix (see GRDFNT above). If mass was 

generated directly (densities in mass units), WTMASS would be 1.0 
and the MO matrix would be a rigid body mass matrix. 

NEW BULK PARAMETERS - PHASE 1 

TINAME label name of INPT, where reduced substructure matrices are out- 

putted for Phase 2. 


MARRC 


■ Converts mass t.tqt #t>i+ ^ 

~ ^eiauj/c incorporated is MaSSC - 

386.4 #/in/sec 5 which converts mass to lbs.^ which is consis- 
tent with the parameter WTMASS = .002588. Therefore, the mat- 
rices MO and MORR will be in consistent units (see GRDFNT above). 
If WTMASS = 1.0, MASSC = 1.0. In order to have MO and MORR 
consistent MASSC should be the reciprocal of WTMASS. 


PHASE 1 ASSUMPTIONS 

1. Any zero- stiffness degrees of freedom and symmetrical or anti- symmetrical 
boundary constraints at the model plane of symmetry are included in the 
Single Fbint Constraint set (SPC). No other degrees of freedom are included 
in this set. 

2. Each substructure should reference the same origin on the GRDFNT parameter 
card and also reference the same basic coordinate system. 

No scalar points should be used in this phase. 
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3. 


4. Interface or junction point degrees of freedom are defined by SUFORT 
cards (r-set). 

5 . The component modes obtained in this phase are with the interface fixed. 
These modes can be plotted. 
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ALTERS INCORPORATED - PHASE 1 (GENERAL FLOW) 



Changes mass at GRID points to weight C w ge"l -CV] * 
MASSC PRINT & CHKPNT 

Default = 386.4 #/in/sec^ 

Obtains from modified module GFWG, basic matrix necessary 
for subsequent built-in checks. Also purges matrices gene- 
rated in ALTER 59, if no MPC’s. 


Performs MPC check. 


Insert 


1 EQUIV D n ,D f /SINGLE 

PURGE CPn^ , D , SPCCK, M , 
f 7 s 7 7 ss 7 

rWCTWV pp^ s T\ T) ciPfPTr M 


w 

ss 


ss 5 

TaT 

? ” 


WS PC/SINGLE 

WSPP 

ss’ 


Equivalences and purges matrices generated in ALTER 65 , if 
no SPC's 


Performs SPC checks 


I EQUIV D , D / OMIT 

PURGE CPf£, D , GOCHK/OMIT 

CHKPNT CPf° D , D , GOCHK 
o’ a’ 

Equivalences and purges matrices generated in ALTER 73, if 
no OMITS 


Performs check on G transformation matrix. 

o 
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- Replaces Rigid Body Modules. The condition that the 
SUFORT card define rigid body supports has been eli- 
minated. The SUPORT card now defines the interface 
d.o.f. This alter also performs a rigid body trans- 
formation check and reduced interface stiffness check. 
Also prints out a rigid body weight matrix to check if 
mass was lost in the reduction process. 

- Obtains free modes fixed at interface 

insert I *“■> » M £ £,,, EED, , CASECC/LAML, M i? 

' OEIGS/C,N,MODES/V,N,NEIGV 

- insert I CHKHW ^ 

OFP LAML, OEIGS, , , , //V,N, CARDNO 

- Generates generalized modal stiffness and mass matrices 
and puts necessary matrices for PHASE 2 on tape INPT. 

- ' Delete CHKHTT statement 

- Replace LAMA by LAML in calling sequence for 3DR2 
module . 


Bl-4 



ALTER 5b - PHASE 1 fDETATT.FT) FTnw') 


r 

V 


ALTER 54 


*NOTE: 


Module GFWG has been modified to output 
matrix EQg. This matrix expresses the 
basic system resultants of the unit g- 
set loads about a chosen origin, defined 
by the input parameter GRDRNT. Actually, 
EQg is equal to the D 1 matrix discussed 
in the NASTRAN PROGRAMMERS MANUAL (sectior 
4 . 29 . 7 ). 



END OF ALTER 5 







ALTER 59 - PHASE 1 (DETAILED FLOW) 
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\ ALTER 65 - PHASE 1 (DETAILED FLOW) . 
ALTER 65 ) =" — - | 



only columns pertaining to symmetric or 
anti- symmetric rigid body motion should 
u6 zero. 


EXTRACT OUT MASS AT SPC’s USING CPn® 


fM | 1 

SS 1 

^[M m ] 

nn 

l 


l 

CHANGE 

MASS M ss TO WEIGHT 

1 1 

w 

M 

1 — 1 

_ II 

» [M ] x MASSC 

SS 

, 



MASSC is 
If none. 


an input BULK parameter, 
default = 386.4 #/in/see'" 


LOSS OF WEIGHT CHECK THRU SPC’s 


[0] = [wspc] = [w ][D ] 
s x 6 SS s 
print 


Col’s represent resultants of weight 
at SPC’s about origin. For symmetric 
or anti-symmetric analyses, only 
col’s for symmetric or anti- symmetric 
resultants should be zero. 


CHKPNT CPn f , D g , D f , SPCCK, M , 


W 


ss’ 


WSPC 


END OF \ 

ALTER 65 / 

. y 
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GENERATE PARTITION VECTOR {CEap, USING VEC INSTRUCTION 


PARTITION RIGID BODY TRANSFORMATION D USING Cl*/ 
— a r 


'V 

"d“~ 
r J 

< — tV 


CHECK TRANSFORMATION MATRIX G 


l 


[0] = [GLCHK] = CG x ][D r ] - [D £ ] 


a x 6 


PRINT 



For symmetric or anti- symmetric analyses, only 
columns pertaining to symmetric or anti- symme- 
tric rigid body motion should be zero. 


REDUCED INTERFACE STIFFNESS CHECK 

[o] 

= [krrck] = [k ][d ] 
r rr r 

r x 6 


^ — PRINT 



CHANGE REDUCED INTERFACE MASS TO WEIGHT 

[W ] = [M ] x MASSC 
rr rr 


RIGID BODY WEIGHT MATRIX CHECK 


[MO ] = [D ] T [W ][D ] 
rr r rr r 

6x6 


PRINT 


MO should be equal to the MO 
rr 

matrix . 

which is printed out from module 
GPWG. For symmetric or anti- sym- 
metric analyses only symmetric or 
anti- symmetric terms should agree. 


J0 

CHKINT CPa r , V> p D r , GLCHK, 

KRRCK, W , MO 

’ rr’ rr 

(end OF ALTER 75, 8^- — 



Bl-10 









PRINT rn 

r^_ n r J ilr.. IT A -i 

LKllJ- L V Lll ti JL ’V 

miNT „ 

[M 1± : = [0p 


[M lr ] = [0,3 [K tr ] 

V — ' ITO/ 


FROM ALTER 75, 84 


TTNTT GENERALIZED FORCES ABO 


[D.f = [EQ 3 = -CD r 3 T [K Ar 3 T [0 je ][K..] 




SUMMATION OF FORCES ABOUT ORIGIN OF MODAL STIFFNESS 
& MASS MATRICES. PRINT OUT ALL MATRICES 
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— I 

COPY ONTO TAPE (INPT) MATRICES NECESSARY TO COUPLE IN PHASE 2 


OUTPUT1 K rp , K^, M rr , M ir , M ±± //C y N, -l/C, N, O/V, Y, TPNAME 
OUTPUT1 D i , D r ,,,//C, N, O/C, N, O/V, Y, TPNAME 


[EXPAND EIGENVECTORS 04 USING C CF&J} j 

- J 2_J_ 




TFNAME is an input BULK Parameter 


-(end of alter 94 ) 

' 


NEW BULK PARAMETERS - PHASE 2 


NOSUB - - - Number of reduced substructures (on tape INPT) to be coupled. 

TPNAME ----- label name of INPT which contains the reduced substructure 


matrices plus column partition vectors for merging. It is 
also the common label name of INP1, INP2, etc., where the 
final substructures system eigenvalues and eigenvectors are 


outputted, which will be used for Phase 3. 

MASSC ----- Same as in Phase 1. 

PHASE 2 ASSUMPTIONS 

1. Interface or junction point degrees of freedom are defined on GRID cards 
with the released DOF in Phase 1 SPC’d out. 

2. Substructure component modes are defined as scalar points. 

3 . Continuity at junction points between substructures are accomplished with 
MPCs. 

4 . SUPORT card in this phase defines the usual rigid body statically deter- 
minate supports. 

5- The Phase 1 tapes are assumed to have been consolidated onto 1 tape by a 
DMAP run, which will be input to this run. This tape also contains the 
column partition vectors necessary for merging. 

6. Free-free modes are obtained in this phase. Plots are not obtained in this 
phase . 
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ALTERS INCORPORATED - PHASE 2 (GENERAL FLOW) 



Deleted (No member input in this phase) 

Contains DMAP loop to combine reduced substructures 


Deleted (No member input in this phase) 


Same as Phase 1 Alters 

NOTE: SUPORT card in Phase 2 defines rigid body 

support s . 


Obtain free-free modes 

| READ K aa , M aa , M rr , G p EED, USET, CASECC/ 
Insert j LAMA, 0a, M., OEIGS/C, N, MODES/V, N, 

( NEIGV 


Revise 


SDR2 (no members defined) 


Insert 


SDR2 CASECC, CSTM, , , EQEXIN, SIL,,,BGPDP, LAMA,, 
0g, , /, jOPHIG, , , /C, N, REIG 
OFP OPHIG, , , , , //V, N, CARDNO 


Extracts out eigenvalues and eigenvectors for each re- 
duced substructure and puts them on separate tapes. 
These tapes are then used in Phase 3. 

Delete error message 


Bl-13 



ALTER 46 - PHASE 2 fDETATTm FLOW) 
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MERGE l St SUBSTRUCTURE REDUCED MATRICES INTO 
PSEUDO- STRUCTURE G- LINEUP, USING COL. PART. 
VECTORS CFGRj & CFGIj 

[ V<= [iff] 


{CPGl] = {CPGIj} 
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^IS \ 
SKLOOP 
\<0 y 


LABEL LD1 


Return to top of loop 
for next substructure 


REPT LOOED, 9 


END OF \ 
ALTER 108, 112) 
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NEW BULK PARAMETER - PHASE 3 

TMAME2 ----- Label name of INPT which contains final substructure system 
eigenvalues and eigenvectors from Phase 2. 

PHASE 3 ASSUMPTIONS 

1. Checkpoint tape from Phase 1 is used in this phase. Also, a tape from 

"HJlo «-»<-»/-* O n /-\V 1 4- ri ■? VT? VI -Pn v> r» 1 nirn4- n«w 

Xi.aa.kjt; (— v-v/uuuxiixx ig, xxixgu. w jr o ocui giguxvc^ uwio, 

2. The final full system eigenvectors are recovered in this phase and can 
be plotted. 


ALTERS INCORPORATED - PHASE 3 (GENERAL FLOW) 



Deleted (done in Phase 1 - on CHKPNT tape ) . 

Deletes part done in Phase 1, and inputs Phase 2 re- 
sults to obtain final a-set eigenvectors. 

Deletes part done in Phase 1. 

Delete error messages. 
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SUMMATION OF ELASTIC FORCES ABOUT REF. PT . 
PRINT 

[0] = [EQFE] = [D ]' 1 [K ][0 ] 

6 x i a aa a 


mNT [SUMM] = [D ] T [M ][0 ] 
c • a aa LK a 
6xi 

For absolute values multiply each column by correspond, 
eigenvalue 

NOTE: Summing up this matrix from all substructures should be 

zero (conservation of momentum). 


ENUALTER 
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